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Superheavy particles are a natural candidate for the dark matter in the universe and our galaxy, 
because they are produced generically during inflation in cosmologically interesting amounts. The 
most attractive model for the origin of superheavy dark matter (SHDM) is gravitational production 
at the end of inflation. The observed cosmological density of dark matter determines the mass of 
the SHDM particle as mx = (a few) x 10 1 GeV, promoting it to a natural candidate for the source 
of the observed ultra-high energy cosmic rays (UHECR). After a review of the theoretical aspects of 
SHDM, we up-date its predictions for UHECR observations: no GZK cutoff, flat energy spectrum 
with dN/dE ~ 1/E 1 ' 9 , photon dominance and galactic anisotropy. We analyze the existing data 
and conclude that SDHM as explanation for the observed UHECRs is at present disfavored but not 
yet excluded. We calculate the anisotropy relevant for future Auger observations that should be the 
conclusive test for this model. Finally, we emphasize that negative results of searches for SHDM 
in UHECR do not disfavor SHDM as a dark matter candidate. Therefore, UHECRs produced 
by SHDM decays and with the signatures as described should be searched for in the future as 
subdominant effect. 

PACS numbers: 12.60.Jv, 95.35. +d, 98.35.Gi 
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I. INTRODUCTION 

Superheavy Dark Matter (SHDM) is an interesting as- 
pect of modern particle physics and cosmology. Being 
first suggested as explanation p], Q for the observation 
of cosmic rays with energy above the so-called Greisen- 
Zatsepin-Kuzmin (GZK) cutoff, it has been later devel- 
oped into the study of a new form of dark matter. 

The concept of dark matter composed by superheavy 
particles with mx > 10 13 GeV (further on we shall call 
them X particles) at first glance seems to be exotic, 
mainly because of two questions: How can SHDM par- 
ticles have a lifetime exceeding the age of the universe, 
and why their abundance should be dominant in the uni- 
verse today? The answers to both questions have been 
known for a long time and are based on quite general 
theoretical concepts. The problem of the stability or the 
quasi-stability of a heavy particle exists also for the light- 
est supersymmetric particle, e.g. the neutralino. Discrete 
gauge symmetries protecting neutralinos from fast decays 
work equally well for other particles. A particular exam- 
ple of a particle with lifetime exceeding the age of the 
universe was given in Ref. 3] . The large abundance of su- 
perheavy relic particles can be provided by gravitational 
production 4] , which works for superheavy particles very 
similar to the production of density fluctuations during 
inflation. 

In contrast to usual thermal relics, like e.g. neutrali- 
nos, SHDM particles are non-thermal relics and have 
been never in chemical equilibrium with radiation. They 
must be produced very early, at the end of inflation. 
Then it is enough to transfer a tiny fraction from the 
energy of radiation to SHDM particles, less than 10 -18 , 
in order to have fix ~ 0.3 now: The energy density of 



non-relativistic X particles diminishes with time as l/o 3 , 
where a(t) is the scaling factor of the universe, while the 
energy density of radiation diminishes as 1/a 4 . When 
normalized at the inflationary epoch, on = 1, a(t) has 
the enormous value ~ 1 X 10 22 now. Not surprisingly, 
this small energy fraction can be transferred to X parti- 
cles by many different mechanisms, such as thermal pro- 
duction at reheating 0, |5| > the non-perturbative regime 
of a broad parametric resonance at preheating ja], and 
production by topological defects 0, Q ■ 

Thermal production of X particles with mx > 
10 13 GeV requires a very high reheating temperature T r h- 
In supersymmetric cosmology, T r h is limited by gravitino 
overproduction as T r h < (10 9 — 10 10 ) GeV. However, the 
gravitino density can be diluted efficiently by entropy 
production during thermal inflation [8(. Thermal infla- 
tion solves the problem of overproduction of particle at 
reheating and allows higher temperatures with efficient 
SHDM production. In Ref. .7], it was shown that the 
maximal temperature after inflation can be much higher 
than Trh- 

X particles are efficiently produced at preheating p|. 
This stage, predecessor of reheating, is caused by oscil- 
lations of the inflaton field after inflation near the mini- 
mum of its potential. Such an oscillating field can non- 
perturbatively (in the regime of a broad parametric reso- 
nance) produce intermediate bosons x, which then decay 
into X particles. The mass of the X particles can be one 
or two orders of magnitude larger than the inflaton mass 
TO^, which should be about 10 13 GeV. 

X particles can be also produced by topological de- 
fects, such as strings or textures. Particle production oc- 
curs at string intersections or in collapsing texture knots. 
X particles can also be produced by hybrid topological 



defects, such as monopoles connected by strings or walls 
bound by strings. The main contribution to the X par- 
ticle density is given by the earliest epochs, soon after 
topological defect formation. Topological defects of the 
energy scale r\ > mx can be formed in phase transitions 
at or slightly before the end of inflation. Efficient produc- 
tion of topological defects is predicted for the preheating 
stage. 

However, the most remarkable creation mechanism for 
SHDM is its gravitational production 9]. Particles are 
created by time-variable gravitational fields during the 
expansion of the universe. For this mechanism the in- 
teraction of X particles with other particles (e.g. with 
the inflaton) is not required, even sterile particles are 
produced. The present abundance f2 s hdm of the X par- 
ticles is mainly determined by its mass mi, while the 
dependence of fi s hdm on the reheating temperature T r h 
is model-dependent. To provide fLhrim = 0.27, needed 
according to WMAP observations E3 , the mass of the 
X particle must be (a few) x 10 13 GeV. 

The stability of X particles can be ensured by discrete 
gauge symmetries. It must be weakly broken, if we want 
long-lived particles with lifetime tx >to, where to is the 
age of the universe. This superweak symmetry breaking 
can occur due to wormhole [lj or instanton effects Q- 
Alternatively, discrete gauge symmetries could be broken 
by higher-dimensional operators 11]. An example of a 
SHDM particle in a semi-realistic particle physics model 
are cryptons, i.e. bound-states from a strongly interacting 
hidden sector of string/M theory [3|, Il2t Il3| . 

What are the prospects to observe SHDM, if X parti- 
cles are absolutely stable? 

This is a pessimistic case for SHDM, because unitarity 
limits severely the AA-annihilation cross section: Since 
the velocity of X particles is very small, t)<l, only the 
s-wave contributes to (ffann") = a+bv +0(v ), resulting 
in an unobservable small cross section, (a ann v) w a < 
l/m x . An interesting and rather exceptional case was 
found in Ref. |l4|, when X particles are cosmologically 
produced in the form of close pairs and form bound states 
due to gauge interactions between them. The lifetime 
of pairs corresponds to the spiral-in time of pairs with 
their subsequent annihilation, as in the case of monopole- 
antimonopole pairs. 

In the framework of gravitational production of 
SHDM, the mass of the X particle is fixed as (a 
few) xlO 13 GeV, and we are left in the SHDM model 
with only one free parameter, the lifetime of the X par- 
ticles tx ■ If one requires that the "AGASA excess" is 
explained by the SHDM model (see below) , then this pa- 
rameter is fixed by the UHECR flux observed by this 
experiment. 

At present, the most interesting manifestation of 
SHDM may be the observations of UHECR beyond the 
GZK cutoff U H El El El El El There are three 
basic signatures of UHECRs from SHDM: 

• SHDM particles as any other DM particles cluster 



gravitationally and accumulate in the halo of our 
galaxy with an overdensity 2.1 x 10 5 . Hence the 
UHECR flux from SHDM has no GZK cutoff in 
energy spectrum [l|. 

• Since in the decays or annihilations of X particles 
pions are more abundant than nuclcons, UHE neu- 
trinos and photons are the dominant component of 
the primary flux [l|. 

• The non-central position of the Sun in the galactic 
halo results in an anisotropic UHECR flux from 
SHDM 0. 

The quantitative predictions for the energy spectra and 
the photon/nucleon (7/-/V) ratio in the decays or annihi- 
lations of X particles required the extension of existing 
QCD calculations for parton cascades from the TeV scale 
up-to the scale mx- The first calculations used the an- 
alytic limiting spectrum approximatio n l2(Jl or extended 
the Monte Carlo simulation HERWIG EH- More recent 
calculations using a SUSY-QCD Monte Carlo JH H2 
and the DGLAP evolution equations[22, l23l l24j pre- 
dict quite accurately the secondary spectra from de- 
cays/annihilation of SHDM particles and agree well with 
each other. Their most important outcome for UHECR 
observations is the flat shape of the energy spectrum. At 
the relevant energies, it can be approximated as dE/E 1 - 9 
while the photon/nucleon ratio is j/N — r 1 i N (x) rs 2 — 3 
[22| . being only weakly dependent on x = lEjnix- 

An anisotropic UHECR flux from SHDM is guaranteed 
by the fact that the distance from the Sun to the outer 
boundary of the DM halo is largest in the direction of 
the galactic center (GC). Numerical simulations of the 
DM distribution show an increase of the DM density to- 
wards the GC as ex r" 1 or oc r" 15 |2Jj. This further 
enhances the expected anisotropy. The relevant calcu- 
lations for this anisotropy were presented in Ref. |17| . 
Comparisons of the calculated anisotropy |18l fl9| with 
existing data of the air-shower arrays on the northern 
hemisphere have revealed no contradiction between data 
and model predictions. By contrast, detectors on the 
southern hemisphere able to observe the GC are much 
more sensitive to this anisotropy. The data of the old 
SUGAR detector located in Australia are only marginally 
consistent with the prediction of the SHDM model J27J • 

At what energy UHECRs from SHDM become the 
dominant component? 

The answer to this question became unambiguous after 
the precise calculation of the spectrum of secondary par- 
ticles produced in the SUSY-QCD cascade, which can 
be approximated as oc E~ 19 . This spectrum is very 
flat and fitting it to the AGASA data J2J| shows that 
it can become the dominant component of the UHECR 
flux only at energies above 8 x 10 19 eV (see Fig.^l. This 
is an important and reliable conclusion about the status 
of UHECR from SHDM. 

The photon dominance is another reliable prediction of 
the SHDM model. Note that this test is relevant mainly 
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FIG. 1: The calculated spectrum of UHECRs from SHDM 
(dotted curves) in comparison with the AGASA data (from 
Ref. |23|). The spectrum from SHDM decays can explain only 
the highest energy events ("AGASA excess"). The dashed 
curve gives the universal spectrum of extragalactic protons 
[2Sj. The sum of these two spectra is shown by full curves 
for two SHDM fluxes with different normalizations. The \ 2 
values are given for the comparison of these curves with the 
experimental data at E > 4 x 10 19 eV. One can see that 
the SHDM flux can dominate only at energy higher than 8 x 
10 19 eV. 



for energies higher than 8 x 10 19 eV, where UHECRs 
from SHDM dominate the flux . Proton and photon in- 
duced showers can be distinguished by the muon com- 
ponent observed at ground level. An analysis of events 
with energies E > 1 X 10 20 eV has been performed for 
the AGASA data [H E2 and very recently in Ref. H2 
for the combined AGASA ^ and Yakutsk data [33|. We 
shall discuss this analysis in Section 2. Here we only note 
that no photon induced showers have been found among 
six AGASA and four Yakutsk events with energies higher 
than 1 x 10 20 eV. One may conclude therefore that the 
SHDM model is not confirmed by this analysis. 

We summarize our introduction emphasizing that at 
present only the "AGASA excess" at E > 1 x 10 20 eV 
motivates the SHDM model as explanation for the ex- 
isting UHECR data. The data of other detectors, e.g. 
Yakutsk or HiRes, are compatible with the GZK cutoff. 
In this case, the UHECR flux from SHDM can be only 
subdominant and is reduced compared to the fit to the 
AGASA data shown in Fig. 1. If the Auger experiment 
confirms the "AGASA excess" , the status of this model 
will be changed. Even in this case, the Auger detector 
has the great potential to confirm or to reject the SHDM 
model for UHECRs by testing its clearest signature, the 
anisotropy towards the galactic center. 

One must clearly emphasize the following. 

The observations of UHECR cannot exclude SHDM 
as explanation for the dark matter in the universe and in 
our galaxy. Assuming that the SHDM are gravitationally 
produced, the mass of the X particles is fixed and the 



only free parameter of the SHDM model is the life-time 
tx- The "AGASA excess" fixes this parameter as tx ~ 
10 20 years. 

/From the HiRes, Fly's Eye and Yakutsk data, that 
are compatible with the GZK cutoff, only a lower bound 
on tx can be derived. Within this lower bound, SHDM 
may still provide subdominant, but observable effects in 
UHECR observations, and some of the showers observed 
at the highest energies could be induced by secondaries 
from X decays. Thus the search for photons coming from 
directions close to the galactic center remains an inter- 
esting task for the Auger detector, even if the "AGASA 
excess" will be not confirmed. 

The discussion above determines the strategy of this 
paper. In Section [H] we obtain more accurately than 
in our previous work [22J the ratio j/N for the SHDM 
model. Note, that the nucleon flux at E < 1 X 10 20 eV is 
given mostly by extragalactic protons and thus our pre- 
diction is determined by the AGASA excess (see Fig. ^) 
and valid only for the AGASA data. We compare our 
prediction with existing analyses of the j/p ratio using 
the AGASA data. The Yakutsk or HiRes data require a 
larger value for tx due to their agreement with the GZK 
cutoff and thus the 7/p ratio derived for the AGASA 
data is in this case an upper limit, not a prediction. In 
Section [3J we calculate the anisotropy of UHECRs from 
SHDM relevant for Auger observations, before we con- 
clude. 



II. RESTRICTIONS FROM PHOTON-INDUCED 
SHOWERS 

In Ref. J23| we have calculated the ratio r 1 / N (x) of 
photons to nucleons in the QCD cascade initiated by the 
decay of a X particle as function of x = 2E/mx- Here we 
shall compute e 1 (E) — (7/tot)^ as the ratio of photon- 
induced showers to the total number of showers at the 
measured energy E using the following set of equations, 
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where the diffuse fluxes J(E) with indices 'shdm' and 
'extr' refer to SHDM and extragalactic fluxes, respec- 
tively. The SHDM fluxes are taken as average over the 
galactic directions observed by AGASA. As extragalactic 
proton flux we use the universal spectrum from Ref. J23 
as shown in Fig. [l] The SHDM spectra jf dm {E) and 
Jp hdm (E) are taken from Ref. [23] (not using the power- 
law approximation). The photon flux J 7 is evaluated at 
the energy XE, where E is the energy determined experi- 
mentally assuming that the primary is a proton. The co- 
efficient A takes into account the differences in the shower 



development between showers initiated by protons and 
by photons. These differences are caused mainly by the 
Landau-Pomeranchuk-Migdal effect . 

Equation J5J normalizes the total flux at E > 1 X 
10 20 eV to the observed "AG AS A excess" (see Fig.[TJ). 
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FIG. 2: Predicted ratio e 7 as function of the observed energy 
for mx = 1 x 10 13 GeV and for different values of the A 
parameter. The solid curve corresponds to A = 1.2, valid 
according to Ref. SHJ for the AG AS A site. The two limits 
in the figure are those obtained in 31] (e = 0.67) and in |32i| 
(e = 0.50) from an analysis of the AGASA data. 

The predicted ratio QJ is shown in Fig. |2 as function 
of the observed energy E for mx = 1 x 10 13 GeV and 
for various values of A. This factor depends both on the 
local geomagnetic field and the detector type and varies 
therefore from experiment to experiment. In Refs. |3(1 
l31| . the effective A was estimated as 1.2 - 1.3 for the 
AGASA site, while Ref. estimated A = 2 for AUGER. 

The analysis of AGASA data at E > 1 x 10 20 eV 
resulted in the following upper limits: e 7 < 0.67 at 
E > 1.25 x 10 20 eV at 95% CL in Ref. [M 
at E > 1.0 X 10 20 eV at 95% CL in RefiB 

Three remarks are in order. 



and e 1 < 0.5 



• The energy calibration of different experiment s by 
the position of the dip in the energy spectrum J35j 
requires a shift of the AGASA energies by the fac- 
tor 0.9 and of the HiRes energies by the factor 1.2. 
These shifts lead to a very good agreement of the 
AGASA and the HiRes data. Such a shift decreases 
further the tension between our prediction of e T and 
the upper limits from Ref. |3ll l32j| . 

• From the Yakutsk and HiRes data we cannot ob- 
tain a prediction for the j/p ratio within the SHDM 
model, instead we obtain only an upper limit. 
These data agree with the GZK cutoff and can re- 
sult only in the lower limit on tx> and thus in the 
upper limit on the UHE photon flux from SHDM. 
On the other hand, J ox tr is fixed independently, and 
thus our calculations give only an upper limit on 



7/p, which is lower than the curves in Fig. [2] ob- 
tained using the AGASA data. 

• Is it possible to use the combined AGASA and 
Yakutsk data for the j/p ratio as constraint on 
the SHDM model for UHECRs? Such a combina- 
tion is possible but requires additional assumptions 
how this combination is performed. One possibil- 
ity is to determine the true UHECR flux averag- 
ing appropriately the AGASA and Yakutsk data. 
Then the SHDM flux would be lower than the one 
shown in Fig. ^ Consequently, the predicted j/p 
ratio would be also reduced, but could now be com- 
pared with the limit 7/p < 0.36 from the com- 
bined AGASA- Yakutsk data |23. Another way to 
include the Yakutsk data in our analysis is to as- 
sume that the lower flux measured by Yakutsk is 
caused by a statistical fluctuation. This assump- 
tion raises the question of the compatibility of the 
measured AGASA and Yakutsk fluxes. The fluxes 
of the two experiments are already incompatible 
at lower energies, 10 19 — 10 20 eV, where the event 
numbers are high. Thus their difference cannot be 
explained simply by statistical fluctuations. Nor- 
malization by the dip J35J results in an energy shift 
of the AGASA and Yakutsk energies by a factor 
0.9 and 0.75, respectively. After this procedure the 
fluxes coincide perfectly. This implies that the en- 
ergies of two Yakutsk events are considerably below 
1 x 10 20 , and hence the upper limit j/p increases 



compared to the one given in Ref. 
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We thus conclude that the obtained upper limits on 
e 7 (see Fig. |3J) do not exclude, but disfavor the SHDM 
model as explanation for the UHECRs. 



III. ANISOTROPY 

The anisotropy in the direction to the GC is the most 
reliable prediction of the SHDM model. Here we shall 
present the detailed calculations of this anisotropy in the 
form convenient for an analysis of the Auger data. 

The angular dependence of the UHECR flux from 
SHDM is given by 
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where r and R are the distances from the Sun and the 
GC, respectively, and fix is the rate of X particle decays 
given by nx(R)/Tx- As distributions of the DM in the 
galactic halo we use the NFW J2£| and Moore et al. |26| 
profiles, 
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with a = 1 and 1.5 for the NFW and Moore et al. profile, 
respectively. We use R s — 45 kpc as obtained in Ref. [18| . 



The distance to the boundary of halo in the direction 9 where Rh = 100 kpc is the size of the DM halo and 
is given by r Q = 8.5 kpc the distance of the Sun to the GC. 



r m ^{0) = r Q cos 9+ yj R 2 h 
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Changing variable r — * R in the integral of Eq. Q we obtain as convenient formula for numerical computations 
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We define the anisotropy A as the ratio of the flux in the direction of the GC within the solid angle Vt to the flux at 
the same energy in the perpendicular direction, 



A( E) = Jf dm (< e, xe) + jf dm (< e, e) + j eKtr (E)n(9) 
y > > [j^(QQo,\E) + j^ a (9o°,E) + j extt (E)]n(e) 
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Here, J s hdm(< &) is the SHDM flux within the angle 9 
relative to the direction to the GC, i.e. within the solid 
angle 0,(6) = 2ir(l — cos 9), and J ex tr is the extragalactic 
proton flux taken at the same energy as Jf| hdm . Explicitly, 



the angular dependence of the fluxes of Eq. 
by 
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The normalized energy-dependent flux J sn dm in Eq. (JJJ 
is obtained by normalization to AGASA excess according 
to Eq. © at E > 1 x 10 20 eV. 

Graphical results of our numerical computations for 
mx = 1 X 10 13 GeV are presented in Fig.Elfor the NFW 
and Moore et al. profiles, while in the Tables [H and [n] 
numerical values of the anisotropy are given for small 9 
and different energies. The anisotropy A(9, E) weakly 
depends on A and is maximized for small 9, but the need 
for a sufficient number of events will make the choice of 
an intermediate value of 9 more suitable. 



IV. CONCLUSIONS 

Superheavy particles are an interesting candidate for 
the dark matter in the universe. They are naturally 
produced in the expanding universe via gravitational 
interactions, when the Hubble parameter H(t) exceeds 
their mass, H(t) > mx- The observed density of DM, 
D, m = 0.27, determines the mass of the particle as 
ra x ~ (a few) x 10 13 GeV. This makes SHDM a nat- 
ural candidate for the production of UHECR. 



TABLE I: Numerical values of A(8, E) in the case of the NFW 
DM density profile. 
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TABLE II: Numerical values of A(0, E) in the case of the 
Moore DM density profile. 
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The SHDM particles (X particles) can be stable (due 
to, e.g., a discrete gauge symmetry) or quasi-stable (due 
to super- weak discrete gauge symmetry breaking). The 
energy spectrum of produced particles is approximately 
a power-law, oc £V 19 . The dominant primary particles 
are neutrinos and photons. The only free parameter of 
the SHDM model as explanation of the "AGASA excess" 
is the lifetime of the X particles, which is determined 
from the UHECR flux as observed by AGASA as t x ~ 
10 20 years. 

SHDM is accumulated in the halo of our galaxy with 
an overdensity of 2.1 x 10 5 and the produced UHECR flux 
thus do not has a GZK cutoff. The production spectrum 
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the dominance of photons and the anisotropy towards the 
galactic center. 

SHDM as a model for UHECR is at present disfavored 
by the following points: 

• The "AGASA excess" is not confirmed by the 
HiRes, Fly's Eye and Yakutsk data, although the 
statistics of all these three experiments is too low 
to state a serious contradiction. The "AGASA ex- 
cess" can be a combined effect of a small system- 
atic error in the energy determination and of the 
low statistics l35t l36l . 
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• Among 17 events with energy E > 1 x 10 20 eV 
detected by all arrays, there is not a single event 
established as a gamma-induced air shower. 

In this paper we calculated the anisotropy, which can be 
reliably tested by the future data of the Auger experi- 
ment. 

The SHDM model has been tuned to explain the 
"AGASA excess" . If this excess is not be confirmed, or 
the predicted anisotropy is not found, it will not exclude 
SHDM as explanation of the DM in the universe. Ex- 
cluding SHDM as explanation of the observed UHECRs 
will put only a lower limit on tx- The production of 
UHE photons and neutrinos by DM in the halo as a sub- 
dominant (for the observed UHECR) effect will remain a 
signature of this model. 



FIG. 3: Anisotropy A(0, E) as defined by Eq. © for the NFW 
density profile (top) and the Moore density profile (bottom) ; 
in both cases A = 2.0 is used appropriate for the surface 
detectors of Auger |34J|. 



oc E 19 can explain only the "AGASA excess" at E > 
8 x 10 19 eV. The two other signatures of this model are 
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